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Abstract 

We construct a long-range Baxter equation encoding anomalous dimensions of composite 
operators in the SL(2\1) sector of M = 4 supersymmetric Yang-Mills theory. The formalism is 
based on the analytical Bethe Ansatz. We compare predictions of the Baxter equations for short 
operators with available multiloop perturbative calculations. 



1 Introduction 



With the discovery of integrable structure^ in QCD [3j HJ E] and maximally supersymmet- 
ric Yang-Mills theory P [7J [H [9] it appears that understanding of strong coupling behavior of 
anomalous dimensions of composite operators is within reach, at least in the latter gauge the- 
ory. At one loop, the dilatation operators in the large— N c limit is identified with a known 
Hamiltonian of a (graded) noncompact Heisenberg magnet with the quantum space in all sites 
corresponding to infinite-dimensional representations of (super)conformal symmetry algebra of 
gauge theory Lagrangians. While at higher loops it is mapped to yet to be determined putative 
long-range spin chain. There are two generic approaches to integrable models, one based on 
(nested) Bethe Ansatz [10] and another relying on the Baxter equation [TT]. While both give 
identical results for models based on representations with highest and/or lowest weight vectors, 
the Baxter framework applies even when the pseudovacuum state in the Hilbert space of the 
chain is absent. For noncompact super-spin chains, the number of eigenstates is infinite for a 
finite length of the spin chain and the analysis of spectra in this approach is advantageous. The 
Bethe Ansatz approach to multiloop dilatation operator in M = 4 super- Yang-Mills theory was 
successfully undertaken in Ref. [T2l [T3] culminating with conjectured long-range Bethe Ansatz 
equations for the PSU(2, 2|4) spin chain [H]. The alternative formulation within the framework 
of the Baxter equation was unavailable due to lack of efficient techniques to work out Baxter 
equations for graded spin chains even with nearest-neighbor interactions. In Ref. [15] we have 
suggested a regular procedure to construction of the Baxter Q— operators for short-range ST(2|1) 
magnet which is straightforwardly generalizable to supergroups of higher rank. The Q— operator 
which determines the energy spectrum of the chain, and thus the one-loop anomalous dimen- 
sions of the dilatation operator in either holomorphic sector of M = 1 super- Yang-Mills theory 
[T6] or in the minimally supersymmetric noncompact sector of the maximally supersymmetric 
gauge theory [17], was shown to obey a second order finite-difference equation — the Baxter or 
TQ-equation — very similar to the one of the bosonic SL{2) chain [U [5]. Recently we have sug- 
gested a long-range generalization of the Baxter equation in the noncompact one-component 
SL(2) sector of M = 4 supersymmetric Yang-Mills theory [18j . A next natural step is to extend 
the formalism to the noncompact graded SX(2|1) subsector [17] of this gauge theory. Since the 
first principle microscopic formalism to build long-range chains is currently unavailable, we reply 
on an effective approach using the analytical Bethe Ansatz |l9j . The latter allows to determine 
the spectrum of transfer matrices from known Bethe equations and as a result to obtain (nested) 
TQ-relations in terms of eigenvalues Q of Q— operators and eliminating the nested Q— functions 
determine the Baxter equations [20] . 



2 SL (2\1) sector and quantum numbers of operators 

The SL(2\ 1) symmetry arises as a reduction of the full superconformal symmetry group SU (2, 2|4) 
of the four-dimensional f\f — 4 theory on the light-cone and operating on the complex scalar field 
X(z) and a single-flavor gaugino ip{z). Gauge theory leads to a particular realization of the 
SX(2|1) algebra on the space of functions of the light-cone chiral superspace Z = (z,6). Both 



Earlier, integrability was found to emerge in reggeon interaction of high-energy scattering amplitudes [HE]- 



fields can be accommodated into a single M = 1 chiral superfield 



${Z)=iX{z) + 6i/>{z), (2.1) 

which arises as a component of the light-cone M = 4 chiral superfield $v=4 PS] 

$Af=4^ 9 A )\ e 2 = e,ez=v = ■■■ + 9 1 9%(Z) , (2.2) 

depending on four superspace Grassmann coordinates 9 , A = 1,2,3,4. Here we have made 
the following identifications X = 0i4 and ■?/> = A 3 with components of the M = 4 fields. In 
the multicolor limit, the sector is spanned by single-trace non-local operators in the light-cone 
superspace 

0{Z 1 , Z 2 ,..., Z L ) = tr{$(Z 1 )$(Z 2 )$(Z 3 ) . . . $(Z L )} . (2.3) 

Expanding these operators in Taylor series with respect to bosonic and fermionic coordinates 
we get conventional Wilson operators of different field contents with an arbitrary number of 
covariant derivatives T> + = d+ acting on X and ip, 

O = tr{^X(0)^ 2 ^(0)^ 3 ^(0) . . . <9^X(0)} . (2.4) 

Since the number of covariant derivatives is not restricted from above, the representations of 
SX(2|1) to which these states belong are necessarily infinite-dimensional. 

The Af = 1 superfield <&(Z) transforms in the infinite-dimensional chiral representation Yj 
of superconformal spin j = I + b = 1 of the SL(2\1) algebra^. Therefore, the L— field operator 
0(Z\, Z 2 , . . . , Z£) belongs to the tensor product (Vj) <glL . The eigenstates of the SX(2|1) spin 
chain belong to this space and they can be classified according to irreducible SX(2|1) represen- 
tations entering the tensor product. The corresponding local Wilson operators are known as 
superconformal operators O. However, to discuss them efficiently, it is convenient to pass from 
the basis of operators to superconformal polynomials \& '. This can be achieved by means of the 
SX(2|1) invariant scalar product, which projects out the superspace operator O to superconformal 
primaries O [15] 

O a = f X{[DZ k ] 3 WjZ^^Z^OiZ,, Z 2 ,..., Z L ) , (2.5) 
J k=i 

where a summarizes quantum numbers of the states and the <SX(2|1)— invariant measure reads 

[VZ]j = EiLLtL f d 2 z [ £ d9{l -zz + 99Y' 1 . (2.6) 



7T 



I*I<1 



Once we find the lowest weight vectors ty a in these representations, the remaining eigenstates 
can be obtained from ty a by applying the raising operators V + and V + of the superalgebra. The 
Bethe Ansatz and Baxter equation give the spectrum of the lowest weights only. Being the lowest 
weights, the eigenstates ty a diagonalize the operators J and J of the Cartan subalgebra and the 
quadratic Casimir operator C2 acting in (Vj)® L 

C 2 ty ct = JJty a , J^ a = (m + L)V a , Jm a = mm a , (2.7) 



2 Throughout this paper, for an exception of a few places, we use notations and conventions of Ref. [15] 
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such that the lowest weights are parametrized by the vector of quantum numbers a = [L, m, m] 
with m and fh being nonnegative integers. These integers define the transformation properties 
of the eigenstates under dilatations L° = |(J + J) and U(l) rotations B = |(J — J), 

* a (\ 2 z, X9) = \ m+m * a ( z , 9) , A" 1 ^ = \ m -^ a (z, 0) , (2-8) 

where (z, 9) = {zk,9k\l < k < L} denotes the coordinates in the light-cone superspace. In other 
words, (m + in)/ 2 defines the scaling dimension of the eigenstates while m — fh defines its U(l) 
charge. 

There are natural restrictions on possible values of the integers m and fh. Examining the 
transformation properties of ^> a under (12.81) . one finds that fh — m > 0. For m — fh = the 
wave function ty a does not depend on 0's and it is a function of z— variables only. Since 
is the lowest weight, it has also to be annihilated by the lowering operators V~ and V~ of the 
algebra. This leads to ty a = 1 or, equivalently, m = fh = 0. Below we will choose it as a 
pseudovacuum state in the nested Bethe Ansatz which does not have any Bethe roots associated 
with it and possesses vanishing energy. It corresponds to the local Wilson operator (up to an 
overall normalization) 

d [Lm =tiX L (0). (2.9) 

The expansion of the single-trace operators (12. 3p in powers of 'odd' variables truncates at 
order L and reads 

L 

O l (Z) = X L (z) + ...+^ L (z)Y[9 k , (2.10) 

fc=l 

where X L (z) = tr \X(z\) . . .X(zl)) and iPl(z) = tr [ip(zi) . . .ip(z L )] are the lowest and highest 
components, respectively. The highest component in the expansion (I2.10p possesses the maximal 
17(1) charge L. However, it is a descendant of the lowest weight vector 

0) = 9 12 9 23 . . . 9 L ^ LX { l\z) ~ V~e x ...9 L x ( l\z) , (2.11) 

which is proportional to a homogeneous polynomial in 0's of degree (L — 1) found from the 
requirement that ^ a should be annihilated by the lowering operators V~ = de k being the 
lowering operator in (Vj)® L . Here Xl \ z ) ^ s a translation invariant function of Zk (with k = 
1, . . . , L) and we use the convention 9jk = 0j — 9k- Thus its i7(l) charge is in — m — L — 1. For 

instance, for L = 2, there is just one lowest state ^ = 9i 2 x < 2\ z ) with X*i\ z ) = X^i 2 ^ ~ z z) 
[2T] . The above lowest state corresponds to the operator possesses the field content tr[X^]. 
The remaining operators of the supermultiplet with different particle content are deduced from 
this one by applying the step-up fermionic operators V + and V + . For x( z i ~ z< i) — z i2i we 
get conventional two-particle conformal operators C[2,n,n-i] = f r [^(^[+]) n 7n°' 1 ' ) (<9[_]/(9[ + ])'?/'] with 
9[±] = d h+ ± <9 2j+ . 

We conclude then that the possible values of integers m and fh are subject to the constraint 
1 < fh — m < L — 1. The sum m + fh, on the other hand, is unrestricted from above since local 
Wilson operators can carry an arbitrary number of covariant derivatives. 

For the low boundary m — m — 1, the eigenstate ty a has a unit U(l) charge and, therefore, 
it is given by a linear combination of 0's with prefactors depending on z— variables only. The 
latter are fixed from the requirement that ty a has to be annihilated by the lowering generators 
yielding 

^ a (z,9)=V- X { l\z), (2.12) 
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with XL \ z ) being yet another translation invariant function of (with k = 1,...,N) and 
V~ = ^2kdkd Zk being the lowering operator in (Vj)® L . Generally, the states with fh — m = 1 
is related to the L— scalar operator with fh derivatives, while fh — m = L — 1 to L— gaugino 
operators with m derivatives as their descendants, i.e., 

%,m,m-i] ~ tr[d™X L (0)} , [LtL+m ^ m] ~ tv[d^ L (0)] , (2.13) 

respectively. The dilatation operator mixes together different components of the sum carrying 
the same number of 0's. A distinguished feature of the two components, X L {z) and 4>l(z), is 
that the dilatation operator acts on them autonomously. For such states, corresponding to the 
so-called maximal helicity operators [3]. 

The case 2 < fh — m < L — 2, is realized for the spin chain of length L > 4. The eigenstate 
ty a carries the U(l) charge equal toro-m and it is given by a homogeneous polynomial in 0's of 
degree fh — m with the coefficient given by z— dependent functions. In distinction to (12.1 ip and 
(I2.12p . these functions are, in general, independent of each other. For L = 4, the lowest weight 
with fh — m = 2 is 

%,m,m-2] = (0120233b + 012024% + 013034%) XlL(z) , (2-14) 

where \ are translation-invariant functions of the coordinates z and we have also introduced 
notations for djk = dj — % 

3 Bethe Ansatz and Baxter equation 

As we pointed out in introduction, the «SX(2|1) integrable spin chains based on R— matrices can 
be solved via either nested Bethe Ansatz pU] or Baxter approach [TT]. The former relies on the 
existence of a pseudovacuum state in the quantum space of the model. It provides a solution 
to the energy spectrum of the model and leads to expressions for the eigenvalues of transfer 
matrices in terms of two sets of Bethe roots. The transfer matrices is the main ingredient of 
the Baxter approach, with the Baxter operators themselves being certain transfer matrices with 
a special-spectral parameter-dependent — dimension of representations in the auxiliary space. 
However, since one lacks a systematic procedure to construct long-range integrable spin chains 
corresponding to gauge theories, one therefore has to resort to techniques which bypass the 
microscopic treatment and rely on general properties of macroscopic systems. The analytical 
Bethe Ansatz method, which is a generalization of the inverse transfer matrix method, was 
developed to determine the spectrum of transfer matrices for closed chains [19] and serves the 
purpose. In this approach, one uses general properties of the i?-matrix, such as analyticity, 
unitarity, crossing symmetry, etc., to derive various properties of the transfer-matrix eigenvalues. 
These properties are used to completely determine the eigenvalues, assuming that they have the 
form of dressed pseudovacuum eigenvalues. We will solely concentrate on the eigenvalues of the 
transfer matrices, not their eigevectors. Thus we assume this approach based on conjectured 
form of nested Bethe Ansatz equations at higher orders of perturbation theory and subsequent 
use of the so-called analytical Bethe Ansatz to find transfer matrices. 

3.1 Short-range magnet 

The SX(2|1) spin chain has in fact three different pseudovacuum states and, as a consequence, 
one can construct three different nested Bethe Ansatz solutions [22l [T5] . We choose however 
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a single nesting corresponding to the first level vacuum built from scalars X and as a second 
(nested) vacuum we use the primary excitation T> + X. The nested Bethe Ansatz equations for 
the superconformal spin j = 1 of the site read 

^ fa — + fa— ra— 1 + (1) fa (1) + 



0,k 



1=1 



0,1 



u 



(1) 



where — u ± ~. A distinguished feature of this noncompact model as compared with conven- 
tional compact spin chains is that the total spin can now take arbitrarily large values and the 
energy spectrum of the model is not restricted from above even for a finite length of the spin 
chain. For specific low values of [L, m, m] the set of transcendental equations (13.11) can be solved 
numerically. These studies show that all first and second level Bethe roots are real. The one-loop 
anomalous dimensions are determined by the first level Bethe roots uo,k only via the formula 
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(0) 



E 



Let us introduce two polynomials parametrized by the first and second level Bethe roots 



(3.2) 



Qo(u) = l[( 



u 



Uo,k) 



Q 



u 



m—ra—1 



n 



u — u 



(1) 

0,fc 



(3.3) 



k=l k=l 

These are eigenvalues of the Baxter operators of the SL(2\1) spin chain. As is well known the 
number of independent Baxter functions depends on the rank of the symmetry group. In the 
present case there are three. However, only two of them are polynomial in the spectral parameter 
u, with the remaining one being a meromorphic function of u. All transfer matrices of the chain 
can be expressed in terms of these polynomials |15j . 

The transfer matrices can be uniquely fixed using the analytical Bethe Ansatz [19] requiring 
that they should be polynomial in u, free from poles at Bethe roots of the first and second level. 
We will presently consider only matrices correspondingto lowest dimension representation in the 
auxiliary space. This consideration immediately yieldaj 



tq{x) 



u 



,Qo{u-i) 
Qo{u) 



Qo(u + i) 
Qo{u) 



1 



Qo ){ 



u 



(3.4) 



It corresponds to the eigenvalues of the transfer matrices with atypical three dimensional repre- 
sentations (1) + in the auxiliary space. One can also construct the transfer matrix with atypical 
representation (1)_ in the auxiliary space and it reads 



T (X) 



[u 



.Qvju + i) 
Qo(u) 



[U 



Qo(u - i) 
Qo(u) 



J Q { o\u- 



(3.5) 



These transfer matrices are polynomials of order L in the spectral parameter u with coefficients 
determined by the conserved charges of the chain, 

L L 

t (u) = (u-) L + qo,k{u-) L ~ k , f (u) = (u + ) L + ^ qo,k{u + ) L - k , (3.6) 



k=2 



k=2 



3 The results of Ref. [15] are reproduced under the following redefinition of the spectral parameter, u — > iu, 
and the chiral Baxter functions, Qz(u) —* Qa(u — |) and Q^(u) 
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and the leading charge determined by the eigenvalues of the quadratic superconformal Casimir 
operator go.2 = <Zo,2 = C2 = rh(m + L). Eliminating the auxiliary Baxter operator from the 
transfer matrices, one finds a second order finite-difference equation for the Baxter function Q 

mm, 

[r (u)f (u) - (u + u~) L } Q (u) = (u + ) L [r (u) - (u~) L ] Q ( u+l ) + ( u ~) L [f (u) - (u + ) L ] Q (u-i). 

(3.7) 

By solving this polynomial equation in the spectral parameter u, one finds the roots of the 
Baxter function in terms of the conserved charges qo,k and go.fc as well as the quantized values 
of the latter. Once found, the Baxter function determines the eigenspectrum of the anomalous 
dimensions at one loop, 

7 {0) = f (lnQo^-lnQoH))' . (3.8) 

A distinguished feature of the lowest X L (z) and highest i/>l(z) components of the superspace 
operator (I2.3P is that the SX(2|1) dilatation operator acts on them autonomously. As we previ- 
ously established, the state X L (z) is a descendant of the SL (2\1) lowest weight vector (12.121) with 
m — m = 1 while 9\ . . . 0li[>l(z) is a descendant of the lowest weight (12. lip with fh — m = L — 1. 
In both cases, the SX(2|1) Hamiltonian effectively reduces to the Hamiltonian of the SL(2) spin 
chain of length L and spins £ — 5 and £ = 1, respectively. 



3.2 Long-range magnet 

We now turn to multi-loop generalization of the SX(2|1) Baxter equation. The starting point 
of our consideration is a generalization of Bethe Ansatz equations for the sector to all orders in 
't Hooft coupling constant g = g YM VNc/ (2tt). Since the Bethe equations were conjectured for 
the full gauge theory in [Hj, we can get the ones for the minimally supersymmetric subsector 
by removing excitations not belonging to the sector in question. One immediately finds that the 
SX(2|1) nested Bethe Ansatz equations take the following form 



n 



x 1 



j I X k 



X, 



(1) + 
k A j 

3=1 X k x j 



n 



4x; x 



4x k xj 



,Xj)-i6(x k ,Xj) 



xT 



n 

i=i 



x i ~ x j 



J 3 
(1) ' 



(3.9) 
(3.10) 



They are written in terms of the renormalized spectral parameter x = \{u + \/u 2 — g 2 ) [12] using 
the convention x ± = xfw 1 * 1 ] and embody a nontrivial magnon-magnon scattering phase factoiQ 9 
which is indispensable to have agreement with string theoretical calculations [221 Elj and four- 
loop calculations of cusp anomalous dimension in maximally supersymmetric gauge theory [25] 
as was recently demonstrated in Ref. [26J, 



6{x±,x 



oo oo 
m=l n=0 



"1) ^2m,2n+l (g) 



Q 2 n+l{Xj) Q,2m[x 



3, 



[X 



±\2m 



[X 



±^2n+l 



(3.11) 



'We slightly changed notations used in Ref. [26] to accommodate 9 for our needs. 
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c, 



[U°°] 



4i 
3i 
2i 



Figure 1: The integration contour in the magnon scattering phase. 
It is expressed in terms of the single-excitation charges Qk(x), 

&(*,) = (3-12) 



/ j' T;'t' ■ PIS) 



and expansion coefficients depending on the coupling constant 

' 9 -Y +n f°° dt J rn(gt)J n (g t) 

.2) J t(e'-l) 
The anomalous dimensions are determined by the first level Bethe roots x k via 

In the zero 't Hooft coupling limit, the long-range equations (13.91) and the anomalous dimensions 
l(9)l ' 9 2 naturally reduce to one-loop Bethe Ansatz (13.1 j) and (13.81) of the previous section. Anal- 
yses based on the SL{2) reduction, i.e., fh — m = 1, of the long-range equations demonstrated 
[271 [28] correct interpolation of cusp anomaly to strong coupling predictions from string theory 
for the energy of dual rotating string configuration on the anti-de Sitter space [291 EQJ ED [32l [33] . 

Now, analogously to the one-loop discussion of the previous section, we introduce two Baxter 
polynomials build from the first and second level Bethe roots 

m m—m—l 

Q(u) = H(u - u k {g)) , Q«(u)= J] (3-15) 

k=l k=l 

which admit perturbative expansion to all orders in coupling constant Up.(g) = u Q ± + 9 2 U\^ + ■ ■ ■ 
and analogously for u k l \g). The transfer matrix can be constructed using the analytical Bethe 
Ansatz [19] requiring pole free structure at Bethe roots of the first and second level. One gets 

r(x) = (a r)V-t«-)2^ (3.16) 

Q(u) 



,A + (x+) Q( M + _ |A_(x+)+iA + (x+) 

Q(u) 



QU{u+) 
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which is a multi-loop generalization of the transfer matrices with (1) + representation in the 
auxiliary space. Here we introduced the convention 



a±0) = a ± (x) - e(x) , 



(3.17) 



for the difference of the trivial dressing factor, represented in terms of the Baxter polynomial 

m ( „i 

n ' 



9' 



3=1 



Axxf 



e % ' 



\<T T (x) 



(3.18) 



with 



a v [x) 



r 1 dt \nQ( V i-gt) ( / 

7-1 7i \ i 



u + gt 



(3.19) 



and the magnon scattering phase factor 



n 

3=1 



e i0(x,xj) _ g6(x) 



(3.20) 



which takes a more involved form 

•i 



e(x) 



9 



dt , Q(- 
m ■ 



2 .0 /•' fi ,^ 



X 



1 



2rci sinh (ttk) 



In 1 



Axx[k + gs] 



Axx[k — gs] 



(3.21) 



The integration contour in the variable k is represented in Fig. [TJ Equation (13.161) reduces to 
the one- loop transfer matrix H3.4f) of the previous section for g = 0. One can also construct the 
transfer matrix with antichiral representation in the auxiliary space, 



f(x) = {x + ) L e 



Q(u) 



(3.22) 



Q(it) 



_ e |A + (x-)+|A_(x-) 



Similarly to Eqs. (13.61) . the transfer matrices (I3.16P and (13.221) are given by the power series but in 
the renormalized spectral parameter with the expansion coefficients determined by the integrals 
of motion which, in turn, admit and infinite series representation in 't Hooft coupling, 



TlX 



X 



E 

fc>i 



Qk(g) 



T(X 



[X 



[X 



Qk(g) 



k>l 



[X n 



(3.23) 



In contrast to one loop, the all-order transfer matrices acquire nontrivial coefficients in front of 
x l-i Thggg can Iqq found explicitly by studying the large-w limit of (13 . 161) and (13. 22H . The first 
subleading asymptotics of the transfer matrices immediately yields q% and q~i, 



n 2 pi rff 



(3.24) 



a 


(Q) 


CI 1 ) 

ryK L ) 


(2) 


(31 


[3,5,3] 


5 


245 
48 


21475 
2304 




To T r 1 

[3, 7, 5] 


133 

IT 


131117 
23040 


1039405829 
99532800 




[4,5,2] 


32±v / 10 
6 


9334±269v / 10 
1728 


48971080±133936lVl0 
4976640 


5(1259661488±33839563v / 10) 905±23 VTd/-(o\ 
286654464 640 ^ \ / 


[4,6,3] 


17 

3 


5005 
864 


658771 
62208 


423834365 443 a(o\ 
17915904 288 



Table 1: Eigenvalues of selected states up to four- loop order in M = 4 super- Yang-Mills theory. 



where 0±(g) is the coefficient of the leading 1/u— asymptotics of the magnon phase (I3.2ip . 

In order to obtain the Baxter equation for the polynomial Q{u), one eliminates the auxil- 
iary Baxter operator Q from the transfer matrices, and finds a second order finite difference 
equation, which we call by analogy with the one-loop case, the long-range Baxter equation 



T(X)T(X) 



f x + x -^A + (x-)+±A + (z-)+|A_(z+)+±A_(z-) 



Q{u) 



(3.25) 



[x 



T(X 



X 



jX e |A+(x-)+|A-(x-) 



+ (x-) L e A ~^ 



Q{u + i) 

Q(u - i) 



The anomalous dimensions to all orders in gauge coupling constant are determined by the equa- 
tion 

-VT~P (lnQ (| - gt) - In Q (-| - gt))' . (3.26) 
i 7r 



Working our perturbative expansion in coupling constant of both sides of Eq. (13.251) and 
anomalous dimension f)3.26p . 

? 2(fe+l) 7 (fc) ; (3>27) 

k>0 



7(5) = ^2g~ v "^'T 



one can explicitly solve for eigenvalues with small m and fh. For the low boundary fh — m — 1, 
we recover the bosonic SL(2) long-range Baxter equation [341 [T8] albeit with a nontrivial magnon 
scattering phase and reproduce all available multi-loop gauge theory calculations [351 S3 EH [38j 
[Ml [25] and predictions based on integrability [121 [131 HQ] . F° r t ne upper boundary fn — m = L — l, 
the solution to the Baxter equation for L = 3 site chains agree with field-theoretical calculations 
of the three-gaugino anomalous dimensions performed in Ref. [3H] and for other values of L with 
the algebraically constructed two-loop dilatation operator in Ref. [4UJ. A few specific eigenvalues 
for up to four loops are displayed in Table [1] with the magnon phase stepping in at order 0(g 8 ). 

The Baxter equation can be immediately used to find the Sudakov behavior of the anomalous 
dimensions for asymptotically large values of the quadratic Casimir C2. As was established in Ref. 
[521 EH], the minimal anomalous dimension is independent of the twist of Wilson operators since 
while go,2 = <?o,2 = C2 takes large values along the trajectory, all other integrals of motion become 
anomalously small such that the spectral curve determining their asymptotics degenerates into 
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one of twist-two operators. An analysis along the lines of Ref. [18] yields equations for the cusp 
anomalous dimension of the lowest trajectory in the spectrum in the large— C2 limit, 

Q(u)Q (u ±i) = e A±{x±) Q(u ± i)Q {u) . (3.28) 

This results immediately implies universality of the cusp anomalous dimension for operators 
with different field content in agreement with independence of the soft gluon radiation of the 
spin of elementary fields it is emitted from. The fine structure of the spectrum determined by 
subleading expansion coefficients depends on the particle content of composite operators and 
deserves a dedicated study. 

4 Outlook 

Presently we have applied the method of the Baxter Q— operator which plays the central role 
in the method of separated variables to minimally supersymmetric subsector of the maximally 
supersymmmetric gauge theory. Within the formalism, the eigenvalue Qo(u) of Q is identified 
with a single-particle wave function which obeys a Schrodinger equation which coincides with 
the Baxter equation. We have constructed a generalization of the Baxter equation for the graded 
SL(2\1) magnet to all-orders in 't Hooft coupling. This putative spin chain arises in the closed 
SX(2|1) subsector of the M = 4 supersymmertric Yang-Mills theory. The TQ-relation was found 
to admit a second order finite-difference form with coefficients determined by "dressed" funda- 
mental atypical transfer matrices. The main advantage of the above construction is that it can 
be straightforwardly generalized to other subsectors of the M = 4 super- Yang-Mills theory based 
on noncompact supergroups of higher rank. This immediately applies to single-trace operators 
with suppressed particle-number changing transitions, with the SX(2|2) subsector of the M = 4 
SYM being its maximal sector. However, above it the length changing effects take place and 
one has to modify the formalism to accommodate them within the Baxter approach. Analyti- 
cal Bethe Ansatz can be naturally used to construct transfer matrices with higher dimensional 
representations in the auxiliary space. Let us point out that the bilinear combination of transfer 
matrices in the left hand side of the Baxter equation is related to the transfer matrix t n /2,& with 
typical 4n— dimensional representation (b,n/2) in the auxiliary space for n = 2, b = [15] . One 
can construct the fusion hierarchy of transfer matrices to determine higher-dimensional transfer 
matrices along the lines of Ref. |41j . A step in this direction has been recently undertaken in 
Ref. 02J. 

Still, to identify the underlying long-range spin chain, one needs the explicit form of the 
Q— operator. Acting on the Wilson operators in the superspace representation (12.31) . the one- 
loop Baxter operator can be realized as an integral operator acting on the positions of superfields 
in superspace, in a close analogy with the dilatation operator which arises as a coefficient in the 
expansion of its kernel in the spectral parameter [15]. The open question remains to find operator 
representation for higher order Baxter functions and ultimately to all orders in coupling. 

This work was supported by the U.S. National Science Foundation under grant no. PHY- 
0456520. 
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